An estimate of QGP viscosity from STAR data on (p mesons 
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In the Israel-Stewart's theory of dissipative hydrodynamics, with a lattice based equation of state, 
where the confinement-deconfinement transition is a cross-over at Tco=196 MeV, we have analysed 
the STAR data on (j) meson production in Au-|-Au collisions at v'^=200 GeV. From a simultaneous 
fit to (j> mesons multiplicity, mean pt and integrated V2, we obtain a phenomenological estimate 
of QGP viscosity, rj/s = 0.15 ± 0.05 ± 0.03, the first error is due to the experimental uncertainty 
in STAR measurements, the second reflects the uncertainties in initial and final conditions of the 
fluid. A host of STAR data, e.g. (j> multiplicity, integrated «2, mean pr, Pt spectra (pr <3 GeV), 
in central Au-|-Au collisions, are consistent with the estimate of viscosity. 

PACS numbers: 47.75.-(-f, 25.75.-q, 25.75.Ld 



I. INTRODUCTION 

Experiments in Au-fAu collisions at RHIC [l], [3, H, 
0], produced convincing evidences that in non-central 
Au-fAu collisions, a hot, dense, strongly interacting, col- 
lective QCD matter is created. Whether the matter 
can be characterized as the lattice QCD [1, 0| predicted 
Quark-Gluon-Plasma (QGP) or not, is still a question of 
debate. A host of experimental data produced in Au-|-Au 
collisions at RHIC, at cm. energy y^=200 GeV, have 
been successfully analysed using ideal hydrodynamics [3] ■ 
Multiplicity, mean pt, pr-spectra, elliptic flow etc. of 
identified particles, are well explained in the ideal hy- 
drodynamic model with QGP as the initial state. Ideal 
hydrodynamics analysis of the RHIC data indicate that 
in central Au-f-Au collisions, at the equilibration time 
Ti w 0.6 fm, central energy density of the QGP fluid is 
Ei w30 GeV/fm-^ 0. 

However, estimate of the initial condition of the fluid 
can not be creditable unless dissipative effects are ac- 
counted for. In recent years, considerable progress has 
been made in numerical implementation of dissipative 
hydrodynamics 1, i, M, MM,M,M d, M M, M 
[Tg . [20I . [2l[. Unlike in ideal fluid evolution, where ini- 
tial and final state entropy remains the same, in viscous 
fluid evolution entropy is generated. Consequently, to 
produce a fixed final state entropy, viscous fluid require 
less initial entropy density (or energy density) than an 
ideal fluid. Amount of entropy generated depend on vis- 
cosity. Viscosity of a strongly interacting QGP is quite 
uncertain. String theory based models (ADS/CFT) giv e 
a lower bound on viscosity of any matter rj/s > I/Att [22| |. 
Perturbative QCD estimates r]/s ^ 1 [2^. At RHIC re- 
gion, Nakamura and Sakai [l^l estimated the viscosity of 
a hot gluon gas as 77/5=0. 1-0.4. In a SU(3) gauge theory, 
Meyer [2^ gave the upper bound 77/s <1.0, and his best 
estimate is 77/3=0. 134(33) at T = 1.165Tc. Given the 
uncertainty (l/47r < 77/s < 1), it is important to obtain 



a phenomenological estimate of viscosity of a strongly in- 
teracting QGP. Only then hydrodynamic analysis for the 
initial condition of the fluid will be reliable. 

For long, strangeness enhancement is considered as a 
signature of QGP formation f26j . Compared to a hadron 
gas, in QGP, strangeness is enhanced due to abundant 
gg —^ ss reactions. Early produced ss, if survive hadro- 
nisation can lead to increased production of strange par- 
ticles compared to pp or pA collisions. Experimental 
data do show strangeness enhancement [23|. Recently, 
STAR collaboration has measured (j) meson production 
in Au-I-Au coUisions at y^=200 GeV Compared 
to pp collisions, in Au-fAu collisions, (j) meson produc- 
tion is enhanced. As noted in [29j . several unique fea- 
tures of (j) mesons (e.g. hidden strange particle, both 
hadronic and leptonic decay, not affected by resonance 
decays, mass and width are not modified in a medium 
[30| etc.) make it an ideal probe to investigate medium 
properties in heavy ion collisions. STAR data [28j appear 
to be consistent with a model based on recombination of 
thermal strange quarks [3l| . As it will be shown below, 
STAR data on mesons are also consistent with hydro- 
dynamic models and enable us to estimate QGP viscosity 
as 77/.S = 0.15 ±0.05 ±0.03. 



II. HYDRODYNAMICAL EQUATIONS 

In Israel-Stewart's theory of 2nd order dissipative hy- 
drodynamics, space-time evolution of a viscous fluid is 
obtained by solving. 



0, 



Dtt^'" = (7r^''- 277V<''u''>) 



(1) 



(2) 
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Eq[T] is the conservation equation for the energy- 
momentum tensor, T'^" = (e -f p)u^u^ — pg^" -f tt^'', 
e, p and u being the energy density, pressure and fiuid 
velocity respectively, tt^^ is the shear stress tensor (we 
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FIG. 1: Minimally viscous fluid {ri/ s— 0.08) simulation for 
temporal evolution of momentum anisotropy in b=7 fm 
Au+Au collision at RHIC. The solid line is the simulation 
result from VISH2+1 [2l|] and the dashed line is the simula- 
tion resuh from AZHYDRO-KOLKATA. 
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FIG. 2: Black circles are lattice simulation Q for entropy 
density. The black line is the parametric representation to 
the lattice simulations. In the inset, the solid and dashed lines 
are p/e in lattice based EOS and in an EOS incorporating 1st 
order transition 7l. 



have neglected bulk viscosity and heat conduction). Eq[2] 
is the relaxation equation for the shear stress tensor tt'^". 
In Eq[21 D — u^df^ is the convective time derivative, 
V<My>.> ^ ^(yt^u" + Vw^) - ^{^.u){g^'•' - m'^m") is a 
symmetric traceless tensor, rj is the shear viscosity and 
is the relaxation time. It may be mentioned that in 
a conformally symmetric fluid relaxation equation can 
contain additional terms [2l| . 

Assuming boost-invariance, Eqs[l] and [5] are solved in 
(r = — z'^, X, y, rjs = ^Inf^f) coordinates, with a 
code '"AZHYDRO-KOLKATA"',^ developed at the Cy- 
clotron Centre, Kolkata. Details of the code can be found 
in [3, [3 ■ In FigUl we have compared the temporal 
evolution of momentum anisotropy Sp = '^r^'^+Tvy y 
a QGP fluid with a calculation of Song and Heinz 21]. 
Initial conditions are similar for both the simulations. 
Within 10% or less, AZHYDRO-KOLKATA simulation 
reproduces Song and Heinz's |2lj result for temporal evo- 
lution of momentum anisotropy Ep. 

III. EQUATION OF STATE 

Equation of state (EOS) is one of the most impor- 
tant inputs of a hydrodynamic model. Through this 
input macroscopic hydrodynamic models make contact 
with the microscopic world. Most of the hydrodynami- 
cal calculations are performed with EOS with a 1st order 
phase transition. Huovinen reported an 'ideal' hy- 
drodynamic simulation with 2nd order phase transition. 
He concluded that the experimental data (e.g. elliptic 
flow of proton or antiproton) are better explained with 
EOS with 1st order phase transition than with EOS with 
2nd order phase transition. However, lattice simulations 
indicate that confinement to deconfinement transition is 
a cross-over at /xb=0 and in the hb range ~30 MeV 
expected at RHIC dlH. In FigH a recent lattice simu- 
lation for the entropy density is shown. The solid line 



in Figl2]is a parameterisation of the entropy density. 

l^ = a+[(3 + jT][l + tanh^-^i (3) 

From the parametric form of the entropy density, pres- 
sure and energy density can be obtained using the ther- 
modynamic relations, 

p{T) = / s[T)ds (4) 

e{T) = Ts-p. (5) 

Generally, in hydrodynamic simulations, hadronic 
phase is approximated by a (non-interacting) resonance 
hadron gas comprising all the resonances below 2-3 GeV. 
As the lattice simulation cover a wide temperature range 
below the cross-over temperature, Tco — 196(3) MeV, 
we choose to use the lattice based EOS (EqI3][5]) both 
in the QGP and in the hadronic phase. The confined 
phase is not a hadronic resonance gas. The trace anomaly 
(e — 3p)/r'* in the temperature range 140-200 MeV is ap- 
proximately 30% less than that of a hadronic resonance 
gas 0. In the inset of Fig[21 the ratio p/e in the lat- 
tice based EOS is compared with p/e in an EOS with 1st 
order phase transition Q, which model the quark phase 
with bag model, and the hadronic phase by the hadronic 
resonance gas. In 1st order EOS, p/e fall sharply near 
the critical temperature. The fall is smoothened out in 
cross-over transition. Lattice based EOS is also softer 
at high temperature but harder at low temperature as 
compared to a first order phase transition. 

IV. INITIAL CONDITIONS 

Boost-invariant solution of ideal hydrodynamics re- 
quire transverse profile of the energy density {e{x,y)) 



3 



% 1.0 

o 

^ 0.9 

Q. 

^ 0.8 
0.7 



Au+Au 

id. hydro. 
- vis. hydro. (iys=0. 08) 
" vis. hydro. (r|/s=0. 16) 




(c) 




50 100 150 200 250 300 350 400 
Npart 



FIG. 3: STAR data on the centrality dependence of (j> meson 
(a) multiplicity, (b) integrated V2 and (c) mean pt are com- 
pared with hydrodynamical simulation of ideal and viscous 
fluid. 



and fluid four velocity {u{x,y)) at the initial time (t.;). 
A freeze-out prescription, e.g. freeze-out temperature Tp 
is also needed. In viscous hydrodynamics, at the initial 
time (ri), additionally required are (independent) shear 
stress tensor {'K^^"{x,y)) components. One also needs to 
specify viscosity coefficient 77 and the relaxation time r^. 
In the present analysis, initial time for the hydrodynamic 
simulation is chosen to be Ti—Q.2 fm. Hydrodynamics 
analysis of photon data in RHIC Au-|-Au collisions sug- 
gest small initial time, ri=0.2 fm [s^]- Even smaller for- 
mation time is suggested in the analysis of J /ij} suppres- 
35| . Fluid is assumed to have zero velocity 



sion at RHIC 

u{x^ y) = at the initial time. The initial energy density 
is assumed to be distributed as Q, 



e(b, X, y) = e40-75iVpQrt(b, x, y) + 0.257Vco;/(b, x, y)], 

(6) 

where b is the impact parameter of the collision and 
Npart and Ncoii are the average participant and collision 
number respectively. The constant Si does not depend on 
collision centrality. Effect of viscosity is enhanced with 
non-zero initial shear stress tensor [15] and we initialise 
it to boost-invariant value. The independent components 
are initialsied as, tt^^ — tt^^ = 2ri/3Ti, tt^^ = 0. For the 
relaxation time , we use the Boltzmann approximation 
= 3rj/2p. Finally, in the cooper- Frye prescription, in- 
variant yield for the mesons is calculated at the freeze- 
out temperature, rF=140 MeV [36j]. 



TABLE I: The fitted values of the initial central energy den- 
sity (si) and temperature (Ti) of the fluid in b=0 Au+Au col- 
lisions, for difi'erent values of viscosity to entropy ratio (ri/s). 
In the last three rows, /N for STAR data on V2 and 

< > are shown. 
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V. RESULTS 

Central energy density (eq) and fluid viscosity (ij) are 
the parameters in our analysis. We assume that viscos- 
ity to entropy ratio t]/s remain constant throughout the 
evolution and simulate Au+Au collisions for i]/s=0, 0.08, 
0.16, 0.2 and 0.25. We fix the central energy density of 
the fluid to reproduce the meson multiplicity in 0-5% 
centrality Au+Au collisions. The corresponding central 
energy density and temperature are listed in tablelH The 
uncertainty in Si is due to the statistical and systematic 
error in STAR measurements. Due to entropy generation 
in viscous dynamics, initial energy density or tempera- 
ture is less in viscous fluid than in ideal fluid. 

In FiglSl in three panels, STAR data [2^ on the cen- 
trality dependence of (j) meson (a) multiplicity (^^), (b) 

integrated V2 and (c) meanpT (< Pt >) are shown. Ideal 
or viscous fluid, initialised to fit (jj meson multiplicity in 
0-5% collisions, reproduces the STAR data on (j> multi- 
plicity in all the centrality ranges of collisions. STAR 
collaboration measured integrated V2 only in 0-5%, 10- 
40% and 40-80% centrality collisions. Fluid dynamics 
overestimate the flow in very peripheral collisions. In 
central or mid central collisions, the integrated V2 is rea- 
sonably well reproduced both in ideal and viscous fluid 
evolution. From FiglSJ;, we also observe that the STAR 
data on < > are not explained unless rj/s <0.25, oth- 
erwise it is overestimated, /N of the combined data 
sets (see tablelj), enable us to estimate QGP viscosity as 
77/s=0.15 ± 0.05. 

The estimate does depend on the assumed initial and 
final conditions, e.g. initial time, initial velocity, freeze- 
out temperature etc. Other conditions remaining un- 
changed, in viscous fluid (77/5=0.16) evolution, < > 
decreases by ^ 10% if initial time increases from Ti=0.2 
fm to 1.0 fm. Simultaneously, increases. Empir- 

ically, < pi^ > linearly increases with 77/ s. Then for 
higher initial time, initial temperature of the fluid will 
be reduced, but 77/s will increase by ~10%. 77/s will 
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FIG. 4: STAR data ,2^] on pt spectra of (j) mesons in differ- 
ent centrality ranges of Au+Au collisions. The dashed black 
and solid blue lines are (j) spectra from evolution of ideal and 
viscous (?7/s=0.15) fluid respectively. 



and we estimate the viscosity to entropy ratio of QGP 
fluid as, r]/s = 0.15 ± 0.05 ± 0.03. The estimate wifl not 
be affected if hard collision contribution in initial energy 
density increases. In 0-5% Au+Au collisions, increasing 
hard collision contribution to initial energy density from 
75%, though multiplicity is reduced, < pl^ > remain 
essentially unchanged. More hard collision contribution 
will increase initial fluid temperature. 

As shown in FigHl <f> mesons px spectra in central and 
mid central Au+Au collisions arc better reproduced in 
viscous (ry/s— 0.15) fluid evolution than in ideal fluid evo- 
lution. We have not shown here, but 77/5=0. 15 is also 
consistent with experimental data on pion, kaon etc. 



VI. SUMMARY AND CONCLUSIONS 



decrease by 2-4% if initially fluid has small transverse 
velocity, Vr = tanh{ar), a=0.02-0.04. With small ini- 
tial transverse velocity < > is increased, multiplic- 
ity remaining unchanged. < pi^ > increase (decrease) 
with decreasing (increasing) Tp. shows an opposite 
trend. If Tp ranges between 130-150 MeV, rj/s will be 
uncertain by -8%. In AZHYDRO-KOLKATA, fluid evo- 
lution is computed with ~5% accuracy leading to ~9% 
uncertainty in rj/s. They all add to the systematic error. 



To summaries, from the STAR data on the centrality 
dependence of (j> mesons multiplicity, mean pp and inte- 
grated V2 , we have obtained a phenomenological estimate 
of QGP viscosity, 77/s = 0.15 + 0.05 + 0.03, the first error 
corresponds to experimental uncertainty in STAR mea- 
surements, the 2nd due to the uncertainty in initial and 
final condition of the fiuid. We conclude that low px < 
3 GeV, (j) meson data in central Au+Au collisions are 
consistent with hydrodynamical evolution of QGP fluid 
with viscosity to entropy ratio rj/s ~ 0.15. 
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